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Abstract 

In this paper we study the relation between long cycles and Bose-Einstein condensation 
in the Infinite- Range Bose-Hubbard Model. We obtain an expression for the cycle den- 
sity involving the partition function for a Bose-Hubbard Hamiltonian with a single-site 
correction. Inspired by the Approximating Hamiltonian method we conjecture a sim- 
plified expression for the short cycle density as a ratio of single-site partition functions. 
In the absence of condensation we prove that this simplification is exact and use it to 
show that in this case the long-cycle density vanishes. In the presence of condensation 
we can justify this simplification when a gauge-symmetry breaking term is introduced 
in the Hamiltonian. Assuming our conjecture is correct, we compare numerically the 
long-cycle density with the condensate and find that though they coexist, in general 
they are not equal. 
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1 Introduction 

Motivated by the path-integral formulation, in 1953 Feynman [1] studied the relation between 
the statistical distribution of particles on permutation cycles and the occurrence of Bose- 
Einstein condensation (BEC). He conjectured that the presence of long cycles is intrinsically 
connected to BEC. Penrose and Onsager developed these arguments and observed that BEC 
should occur when the fraction of the total number of particles belonging to long cycles is 
strictly positive [2]. These concepts, which are now generally accepted, were made mathe- 
matically precise by Siito [3] who also proved the equivalence between the Bose-condensate 
density and the density of the number of particles on long cycles in the case of the free 
and mean- field Bose gas (see also Ueltschi [4]). Subsequently it was shown that this relation 
holds for the perturbed mean-field model of a Bose gas [5]. In our previous paper the validity 
of this hypothesis was tested in another model of a Bose gas, the Infinite-Range-Hopping 
Bose-Hubbard Model with Hard Cores [6]. There it was shown that while the existence of 
non-zero long cycle density and BEC coincide, these densities were not necessarily equal. 

This paper is the sequel to [6], where now the hard-core interaction is replaced with a finite 
on-site repulsion to discourage but not forbid multiple particle occupation of individual sites. 
The thermodynamics of this model have been studied by Bru and Dorlas [7]. They show 
that the phase diagram for this model is much more complicated than in the hard-core case 
as for low enough temperatures there are several critical values of chemical potential which 
correspond to intervals of BEC (see Fig 2). As in [6], we use standard properties of the 
decomposition of permutations into cycles, to convert the grand-canonical sum into a sum 
on cycle lengths. This makes it possible to decompose the total density p = pshort + Piong into 
the density of particles belonging to cycles of finite length (pshort) and to infinitely long cycles 
(piong) in the thermodynamic limit. We consider the relationship between Bose-condensation 
and long cycles in this model. 

The model is considered in the grand-canonical ensemble. In terms of the random walk 
representation, the particles hop from one site to another with probability depending on the 
occupation number of the destination site - the more particles on the site, the less likely 
another particle will hop there. Following [6] we write the cycle density for the number of 
particles on a cycle of length q in terms of a partition function of q distinguishable particles 
interacting with the Boson system. Again as in [6] we exploit the fact that the hopping 
between the q distinguishable particles can be neglected, to obtain an expression for the 
cycle density in terms of the ratio of the partition function for a Bose-Hubbard Hamiltonian 
with a single-site correction and it without. Inspired by the Approximating Hamiltonian 
method [8] we conjecture a simplified expression for the short cycle density Pshort as a ratio 
of single-site partition functions. We prove that this simplification is exact in the absence 
of condensation and implies a zero long-cycle density in this case. Unfortunately we are not 
able to prove this simplifying conjecture when BEC occurs, however we can go some way 
towards justifying it by introducing a gauge-symmetry breaking term in the Hamiltonians. 
Assuming our conjecture is correct, we perform some simple numerical techniques to compare 
the long-cycle density with the condensate. We find (as in [6]) that though they coexist, in 
general they are not equal. 

Before describing the layout of this paper, it is worth noting that BEC may be classified into 
three types (see [9] and [10]): type I/II when a finite/infinite number of one-particle quantum 
states are macroscopically occupied (resp.), and type HI when no states are macroscopically 
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occupied. The relation between the size of long cycles and these condensate types for the 
free Bose gas is considered in [11]. 

The paper is structured as follows. In Section 2 we first describe the model and recall its 
thermodynamic properties as stated by Bru and Dorlas [7]. In Section 3, by applying the 
general framework for cycle statistics described in [5] (following [12]), we form an expression 
for the density of cycles of length q by isolating q distinguishable particles from the boson 
field and show that we can neglect the hopping of these q particles in the thermodynamic 
limit. Section 4 deals with the above conjecture, proving its correctness in the absence 
of condensation and proving an equivalent result with the addition of a gauge-symmetry 
breaking term to the Hamiltonian which is correct in the absence and presence of BEC. 
Section 5 proves that in the absence of BEC the long cycle density is zero. 



2 The Model and Results 

The Bose-Hubbard Hamiltonian is given by 

Hf' = J Yl ia:-al){a,-ay) + Xj2Mn.-l) (2.1) 

x,yGAv ■.\x—y\ = l xeAy 

where is a lattice of V sites, a* and are the Bose creation and annihilation operators 
satisfying the usual commutation relations [a*, tty] = 6x,y and rix = a^a^. The first term with 
J > is the kinetic energy operator and the second term with A > describes a repulsive 
interaction, as it discourages the presence of more than one particle at each site. This model 
was originally introduced by Fisher et al. [13]. 

The infinite-range hopping model is given by the Hamiltonian 

^^ = ^Y. -%) + ^Y. - (2.2) 

This is in fact a mean-field version of (2.1) but in terms of the kinetic energy rather than the 
interaction. In particular, as with all mean-field models, the lattice structure is irrelevant 
and there is no dependence on dimensionality, so we can take = {1, 2, 3, ... , V}. The 
non-zero temperature properties of this model have been studied by Bru and Dorlas [7] and 
by Adams and Dorlas [14]. Also Dorlas, Pastur and Zagrebnov [15] considered the model in 
the presence of an additional random potential. 

Bru and Dorlas applied the "Approximating Hamiltonian" method (see [8, 16]) to the 
Infinite-Range-Hopping Bose-Hubbard Model. In this method one performs the following 
substitution for the Laplacian term of the Hamiltonian: 

1 ^ ^ 

x,3/=l x=l 

(some r G C) to obtain the approximating Hamiltonian: 

V V V 

Hf''{r) = - $^(™. + ral) + V\r\' + A ^n.(n. - 1). (2.3) 

x=l x=l x=l 
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Introduce a gauge breaking source z/ G C in both Hamiltonians (2.2) and (2.3), by setting 
Hy{u) := Hy - ELii^a, + and //^^^(r,//) := H^^^{r) - Y^Lii^^, + K)- Then 
for all /X and z/, one finds that the pressures for these Hamiltonians are equivalent in the 
thermodynamic limit, i.e. for large V one obtains the estimate 

< Pv[Hy{u)] - supp,[i/^PP(r, u)] < 0{V-'/') 

reC 

where for a Hamiltonian H, Pv[H] := py[H]{P, /i) denotes the corresponding grand-canonical 
pressure. Henceforth the /3 and /i dependencies are assumed unless explicitly given. 

With this technique, Bru and Dorlas managed to obtain the limiting pressure and showed 
that in some regimes Bose- Einstein condensation occurs. They proved the following result: 



Theorem 2.1 The pressure in the thermodynamic limit for the Infinite- Range- Hopping 
Bose-Hubbard Model, p{/3,fj,) := limy^aoPviHv], is given by 

p(/3,/i) = sup ( - + ilntrace^^(c)e-^'^«| (2.4) 
r>o y P ) 

where 

h{r) := (1 — fi)n + Xn{n — 1) — r(a + a*) 

is a single site Hamiltonian with creation and annihilation operators a* and a, and with 
number operator n = a*a. Note that it is sufficient to take the supremum over the set of 
non-negative real numbers. 



The Euler- Lagrange equation for the variational principle is 



I ^> tracejr(-c)(a + a*)e ^^^'''^ 



Moreover the density of the condensate is exactly given by 

V 



x,y=l 



where is the largest solution of {2.5). 



Equation (2.5) can have at most two solutions. Clearly r = is always a solution. When /3 
is large enough, for certain values of a second non-zero solution may appear (see Fig 1). 
So := unless a second solution r > exists, in which case : = r. 

The properties of this model were then obtained numerically by finding this maximal solution 
of the Euler-Lagrange equation and then evaluating the pressure using (2.4). As may be 
seen from Fig 2, for sufficiently large /3, there may exist several critical values of which 
correspond to intervals of = and > 0. 

In addition Bru and Dorlas showed that Theorem 2.1 holds in the presence of the gauge- 
symmetry breaking term. In that case, the corresponding Euler-Lagrange equation has a 
unique non-zero solution r^{u). 



Figure 1: Comparison of 2r with {a + a*)h{r) with /3 = 4, A = 5 for the cases fi = 1.5 
(condensation) and /i = 5 (no condensation). 



In this paper we shall analyse the cycle statistics of this model. Hy is a grand-canonical 
Hamiltonian given by (2.2) acting upon the bosonic Fock space T+iHv). Let i/^"^ be the 
restriction of Hy to the n particle space 'H^\ and its corresponding symmetrised subspace 
'H^^^. Then the grand-canonical partition function for this model may be written as 

oo 

(n) 



n=0 



E^, = trace ^^(^,) [e-'^^^-"'^^-)] = trace 

oo 



n=0 7r65n 



where Sn is the set of all permutations of n items, and the unitary representation of a 
permutation tt ^ Sn- 

There is a natural probability measure (see [5]) on the set of all permutations IJ^o (taking 
5*0 = {1}) defined as 



= ^ V ^trace^(„) f/.e-^(<'-^") 



where X is the indicator function, ensuring n ^ Sn for some n. 

From the random walk formulation (see for example [17]) one can see that the kernel of 
Q-p^v is positive and therefore the righthand side of this expression is positive. 

Each permutation n E Sn can be decomposed uniquely into a number of cyclic permutations 
of lengths gi, g2, • • • , Qj with j < n and gi + g2 H — ' + Q'j = For ^ {I5 2, . . . , n}, let Nq^n) 
be the random variable corresponding to the number of cycles of length q in vr. Then the 
expectation of the number of g-cycles in the grand canonical ensemble is 

00 
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Figure 2: Plot of the condensate density r^{= ) versus fi, for /3 = 4, A = 5. 
and the average density of particles in g-cycles is 

This brings us then to the following definition. 

Definition 1 The expected density of particles on cycles of finite length is given by 



Pshort 



lim lim y'c(^(g) 

1— >oo V^oo ' 



(2.6) 



9=1 



and the expected density of particles on cycles of infinite length is given by 



r long 



lim lim V <(g). 

1— >00 V^OG ' 



(2.7) 



g=Q+l 



Clearly = p^^^^^^ + p^ 



ong 



For brevity, denote c^{q) = limy^oo Cv{q)- It is clearly easier to deal with p^^ort since we can 
take the thermodynamic limit inside the sum over q to get 



oo 



Pshort = 2^ 
9=1 
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For the free Bose gas, the mean-field and the perturbed mean-field Bose gas, it has been 
shown that = the condensate density. However in the case of the Hard-Core 

Infinite- Range- Hopping Bose-Hubbard model (see [6]), a different conclusion was obtained: 
that pj^j^g > if and only if p'^ > 0, but that in the presence of condensation these quantities 
were not necessarily equal. We wish to argue that this is also the case for the chosen model. 

We now state our results. 

Theorem 2.2 The density of cycles of finite length in the IRH Bose-Hubbard model may be 
expressed as 



-(S{q-i,)q jj^^ trace ^^(^y) exp { - /3(2Agni + Hy - pNy)} 

V-inr, fvQno ^ , s avr, ! — R( TJ — 1 1 AT \\ 



Pshort " v^oo trace t+chv) exp { - (3{Hy - pNy) } 

where rii is an operator which counts the number of bosons on the site labelled 1. 

If one were to substitute Hy for the approximating Hamiltonian Hy^^{r^) (where again 
is the maximal solution of the Euler- Lagrange equation (2.5)) into the right hand side of 
this expression, one would obtain: 

a trace ^^(c)e-^^^('^^) 

Pshort trace ^^(oe-z^'^oCrM)- ^^'^^ 

where 

hq{r) := (1 — p){n + q) + \{n + q){n + q — 1) — r{a + a*) 

is another single-site Hamiltonian (note that hQ{r) = h{r)). This leads to the conjecture 
that (2.8) gives the correct expression for pshort- 

Moreover the fact that a state corresponding to Hy in the thermodynamic limit may be 
shown to be a convex combination of one-site product states of the form 

trace j-_^(C)e-'^'*«('>M 
^ ' trace j-^(c)e-^''o(^A.) 

supports this conjecture. We prove the conjecture for those values of p such that = 0, 
but unfortunately are unable to do so when > 0. However we can prove a slightly weaker 
result with the addition of a gauge-symmetry breaking term. 

Let Cv(g, v) be the density of particles on cycles of length q for the gauge-symmetry broken 
Hamiltonian Hyiv). 



Theorem 2.3 For p such that = 0, we have 

c\q) 



trace jr^(c)e ^''"("^ 



trace j-+(c)e-^''o(°) ' 

More generally for any p G M, for a fixed u > there exists a sequence Uy u as V ^ oo, 
independent of q such that 

, trace ^^(ae-^['^''('^''(''))-''("+"*)] , , 

lim < g, Uy) = r r ^ 2.9 

v^oo trace J- (c)e-/3['»o('-M(f))-f(a+a*)] ^ j 
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Note that this theorem imphes (2.8) for = 0. This allows us to show that in the absence 
of condensation all particles are on short cycles. 



Theorem 2.4 In the absence of condensation, i.e. for those /i such that 
of particles on short cycles equals the total density of the system, that is: 



0, the density 



Pshort 



Pi 



ong 



0. 



Moreover considering equation (2.9), if the thermodynamic limit and the limit removing the 
gauge-breaking source are interchangeable then (2.8) follows for any r^. 

Assuming the conjecture is correct, we apply some simple numerical techniques to (2.8) 
in order to compare long cycles with the Bose-Einstein condensate. As may be seen from 
Figures 3 and 4 the calculations certainly agree with Theorem (2.4), i.e. that the absence 
of condensation implies the lack of long cycles and visa versa. However more importantly 
they also indicate that while the presence of condensation coincides with the existence of 
long cycles, their respective densities are not necessarily equal. In fact, one can see that 
the long cycle density may be greater than or less than the condensate density for differing 
parameters. 



P=1 , >.=5 




Figure 3: Comparison of p^^^^ with p'^ 
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Figure 4: Comparison of p^^^^ with for various values of (3 and A. 
3 Proof of Theorem 2.2 

Before proceeding to the study the cycle statistics for this model we need to define the n- 
particle Hamiltonian in more detail. The Hilbert space for a single particle on a lattice of V 
sites is "Hv- •= and on it we define the operator 

K = I- Py 

where Py is the orthogonal projection onto the unit vector 

1 ^ 1 
= ^ V = ^(1, 1, . . . , 1) e 

with {ex}^^i the usual orthonormal basis for l-Ly. hy is the orthogonal projection onto 
the subspace orthogonal to gy. For an operator A on 1-Ly^ we define A^"^^ on 7/^^ : = 
® "Hv ® ■ ■ ■ ® ^v, by 

v ' 

n times 

A^"-^ = A®I®...®I + I®A®...®I + ... + I®I®...®A. 

Let J^{'Hv) = ^^^^ denote the unsymmetrised Fock space of l-Ly and define dT[A) on 

T{Hv) as 

oo 

dT{A) = 

n=0 
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where 7/^^°^ := C and H^^^^ ■.= 'Hv. With this notation we can write the free Hamiltonian 
acting on J^{l-Ly) as: 

H'^'^' = dT{hy). 

This represents a collection of particles on a lattice of V sites which hop freely from site to 
site with no inter-particle or external interactions. The hopping action is reflected by the 
Pv operator. 



For bosons we have to consider the symmetric subspace of J^{'Hv)- The symmetrisation 
projection a" on "H^^ is defined by 

< = 1 f/. (3.1) 

where f/^ : i— )■ Ti^P^ is the unitary representation of the permutation group Sn on 
defined by 

U^i(pl O 02 ® ■ ■ ■ ® 0n) = 07r(l) ® 07r(2) ® " " " ® 0,r(n), (pj ^ "Hy, j = 1, . . . , Ti; TT G S'„. 

The symmetric n-particle subspace is "H^^ := (t'^'H^\ allowing us to define the symmetrised 
Fock space as J^+CHv) ■= 0^=0"^^"+- 

The operator which counts the number of particles at site x, n^, is defined by Ux = dT^P^) := 
dr{\ex){ex\). Then the total number operator is A''^ = X]I=i'^a;- 

Let us define the Hamiltonian Hy on J-'{'Hv) by 

V 

Hy = dT{hy) + Xj2nAn.-l). (3.2) 

X=l 

This Hamiltonian restricted to J-'+{Hv) is in fact the Bose- Hubbard Hamiltonian. 

The proof of Theorem 2.2 is in three steps. First we obtain a convenient expression for 
d^{q), the density of particles on a cycle of length q (Lemma 3.1). This involves the par- 
tition function of q distinguishable particles interacting with the boson system through the 
Hamiltonian (3.2). Then we construct a modified cycle density, denoted Cv{q), which ne- 
glects the hopping of the q distinguishable particles and show that these cycle densities are 
equivalent in the limit (Lemma 3.2). Finally we simplify Cv{q) (Lemma 3.3). 

We shall denote the unitary representation of a g-cycle by Uq : — > "Hv"*, that is 

Ug{(f)i, (g)---(g) = (g) . . . (f)i^ (g) (pi^. 

Denote the identity operator on l-iy^ by J*^'^\ and upon J^^{'Hv) by I . When there is no 
ambiguity we shall simply write Ug for Ug^I : J^+{'Hv) — > J^+{'Hv)- Note that [Ug, cr"] = 0. 

We state three lemmas without proof in the course of the argument of Theorem 2.2 and 
prove them shortly afterwards. 

Lemma 3.1 The density of particles on cycles of length q is 

<{q) = ^trace^,,, [[/,e-^(^-'^^-)] 
where '^y is the grand- canonical partition function for Hy and 'Hq^y '■= 'Hy '' ® J-'+{'Hy). 
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This lemma uses cycle statistics to split the symmetric Fock space J^+{'Hv) into the tensor 
product of two spaces, an unsymmetrised g-particle space "H^ ^ and a symmetrised Fock space 
J^+iUv). Write 

A^") = ® I and dT\A) = J^*^) ® dr{A) 

for any operator A on T-Ly. In this fashion, the number operators apphed to 'Hg^v are defined 
as 

= ® I and = /^^^ ® dT{P^). 

Then we may define Hy^ on "H^^v by 

V 

# = /il'^ + dT'ihy) + X ^(n, + N,){n, + N, - 1). 

Define a modified Hamiltonian which neglects the hopping of the q distinguishable particles 
as follows: 

V 

#(^) = + dT'ihy) + X Y^i^X + iVx)(n. + iV. - 1) 

a;=l 

SO that i/v'^'' = Hv"^ — Py \ and define the corresponding cycle density (henceforth called the 
"modified cycle density") by 



Then we have the estimate: 
Lemma 3.2 



l<(g)-?(;(g)|< 4^ ^. 



This implies that in the thermodynamic limit, we are able to disregard the hopping of the 
g-unsymmetrised particles in the cycle density. 

The modified cycle density can be re-expressed as: 
Lemma 3.3 

trace T+inv) exp { - P{Hv - ^Ny) ] 
where rii is an operator which counts the number of bosons on the site labelled 1 of the lattice. 

Combining the above information, we deduce that 

= e-/3(.-M)a liin trace ^^(^^) exp { - (3{2Xqm + Hy - /iiV,.)} 
v^oo trace T+^-Uy) exp { - l3{Hy - fiNy) } 
and this completes the proof of the Theorem. 

Now we shall prove the lemmas. 



11 



3.1 Proof of Lemma 3.1 



The canonical expectation of the number of g-cycles may be found to be 



trace , 



(n) 



(?7g® )e-^^v' 



by following the proof of Proposition 3.1 of the preceding paper [6] and omitting all the 
hard-core projections. 

Then going to the grand-canonical ensemble we obtain: 

e^^'''Zp{n,V)E'^{Ng) 



n=q 

oo 

E 

n=q 
1 



^ trace 



q 



trace 



Hence 



as desired. 



V 



1 

trace ^{q) 



-/3iHv~^lNv)] 



□ 



3.2 Proof of Lemma 3.2 



The following technique has been employed to prove a similar result in [6]. However in this 
case there are several important differences and therefore we give the proof in full. 

To prove Lemma 3.2 we have to obtain an upper bound for 



trace ^ 



<1,V 



JJ ^-li{Hv-t^Nv) 



trace 



In order to do this we first shall introduce some notation. Let {0fc}^o orthonormal 
basis for J^^iTiy)- 



Let Ay^ be the set of ordered g-tuples of (not necessarily distinct) indices of and for 

,^2, • • ■,iq) 



i = Hi, 12, . . . ,i„) e Av^ let 



\i) = \ii,i2, . . . , iq) = ei^ (g) (g) ■ ■ ■ (g) e^^. 
Then {|i) |i G aI^^} is an orthonormal basis for "H^^ 

A basis for l-iq^v may therefore be formed by taking the tensor product of the bases of "H 
and J^+{Uv), so the set {|i) (g) 0/, | = 1, 2, . . . ; i G A^^^} is an orthonormal basis for Hn y. 
For brevity we shall write 

\i-k) = \i)®<Pk- (3.3) 



(9) 
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For simplicity, denote P := Py^ and H := Hy^ — fiNy We expand 

tracer,,, [U,e-f^^"^->^^^^ = trace [U,e-^^^-^^ 
in a Dyson series in powers of P. If m > 1, the m^^ term of this series is 



/"I rsi rsm-1 

= Z?™' / dsi / ds2 ■ ■ ■ dsm tracc 
Jo Jo Jo 



(3.4) 



Let Pr = I ® 



r*''place 



/, SO that P = Pr. Then 



= Z?'" 5^ 5^ ■ ■ ■ 5^ ^n^(rl, r2, . . . , r„,) 



ri=lr2=l r-m=l 



where 



1 pSi pSm-l 

, , 2, . . . , , my — / c^si ds2. . . / rfsm trace 
-'0 Jo 



jjH{l-si)p -pH{si-S2)p 



(3.5) 



In terms of (3.3), the basis of "Hg^v, we may write 



dsi ds2... / XI ■ ■ ■ 5Z 

■ ■ ■ (i™^^ A;™-i|e-^^(''"-^^"'")p,„|i™; A;"^|e-^^''"t/g|i°; (3.6) 

where it is understood that the i summations are over A^^"*, the set of ordered g-tuples (not 
necessarily distinct) of A^^, and the k summations are over the bases for J~j^{l~iY^. 



Notice that we may express 



e-^^^|i;A;)=e-^'?(i-'^)^|i;e-^^'^|A;) 



where 



V 



H' = Hy- ^iNy + A Y^iKiNl - 1) + 2A^>.) 



X = l 



and A^^ = Yl'j=i counts the number of particles at site x which are in g-space. Also, for 
any fixed r: 



1 ^ 

Pr|i; k) = — . . . ,ir, j, ■ ■ ■ , iq, k) 



(3.7) 



where the hat symbol implies that the term is removed from the sequence. 
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It is convenient to define the operation [r, x](i) which inserts the value of x in the r^^ position 
of i instead of v- So for example taking the ordered triplet i = (5,4,1), then [2,8](i) = 
(5, 8, 1). We shall denote the composition of these operators as [r^, x^; . . . ; r2, X2; ri,xi] : = 
[rk,Xk] o ■ ■ ■ o [r2,X2] o [ri,xi]. 



So (3.7) may be written as 

^r|i;fc) = :^X^|[r,j](i);fc). 



Using these facts, a single inner product term of (3.6) may be expressed as 

(i;fc|e-^^^P,|j;A;') =6"^'^'"''^' {k\e-^"''\k') (i|P,|j) 

-/3<?(i-m)s . J" 

-(A;|e-^^'1A:')5^(i|ji,...j„m,...,j,) 

m=l 
V 



V 








V 









{k\e-^''''\k')Y,^h,n■■■^^r,jA.,m...S.. 



m=l 



V 



Now if we sum over j 



^(i; k\e-^^^P,\y, k'){y, k'\ = (A;|e-^^'lfc') ^(^i, . . .X,Jr, • • • k'\ 

j ir = l 



ir = l 

Performing two summations for fixed ri and r2 we get: 

„-/3g(l-At)s .0 
-f}q{l-ti)(s+t) .0 ,, , ,,=0. 



y2 



?1 =1 ^2 ^1 
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Thus (3.5) looks like 
Xyn{ri,r2, ...,rm) 

„-/3(?(1-m) /•! /-SI 



Jo 



V V 



dsi / (is2 ... I ds 



^0 
V V 



m — 1 







E E E ■■■ E ([^-C; ... ;r2,^.^,;n,^^J(i°)|f/,i°) 



:0 ,-1 —1 ,-2 _-| jm =1 



trace 



From the Holder inequality (see for example Manjegani [18]), for non- negative trace class 
operators Ai, A2, . . . , Am+i we have the inequality 



m+l 



|trace (^41^42 . . . Am+i) \ < trace |v4iy42 . . . A^+i \ < Y\_ (t^^-ce A^''^ 

k=l 



j_ 

Pk 



where ^^^^ - = 1, > 1. 

'1 . . . ; Pm 



Set pi = p2 = ^ 
trace 



-^3 — ,Pm+i = —. Taking the modulus of the above 



trace 



-I3H' 



1-si 



< trace jr^^nv 

■ .. trace jr^(^Hv) 



trace 



Since the trace is independent of the V — q sites {i", [ri, zjj(i°), . . . , [r^, 
; r2, 2^2! '^i' '^ri](i°)}5 the product of all the trace terms above is equal to 

trace j-^CH^) 
with \ = {V - q + l,V - q + 2,...,V}. 

This is independent of the i° and i summations, so we need only consider 



V V 



V 



E E ■■■ E E([^-'C; ••• ;-2,v^,;ri,z^j(i°)|t/,i°). 



(3. 



il =1 =1 
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Fix the values of ri, . . . , and zj^, i^^, . . . , zj?^. We intend to show that 

5^([r„,z™; ... ;r2,^?,;n,^JJ(i°)|f//) = l. 



If {ri,r2, . . .,rm} 7^ {1,2, ...,q}, then |[r,„,2™ ; . . . ; r2, 2^^; n, is of the form 

where {ni,n2, ■ ■ ■} is a non-empty ordered set of distinct integers between and q. This 
vector is clearly orthogonal to Ugi^ except for the single choice of 



^ |j2! • • • ) Jni — 1) Jni! • • • ^ inn jr. 



For the case {ri, rg, . . . , r,„} = {1,2, notice that \[rm,iZ,y ■■■ ; ^^2, ^^3! '"i' 

independent of i*^ so we may take it to be 



where s" = (1, 2, 3, ... , q). For each choice of z^^, i^^, . . . , there exists only one possible 



i° G Ai^^ such that 



So we may conclude that 

V V 



V 



(3.9) 



:0 ,-1 —1 ,-2 —1 jm 
I l^^ —1 '■m 



and by using this, we see that the modulus of (3.5) may bounded above by 



|X„(r„r2,...,rJ| <trace^^(^,) le'^^'l e-^-^^i"^) 



V V 



V 



:0 ,-1 —1 ,-2 _i im —I 



trace 



-(3H'- 



ml 



which is independent of ri, r2, . . . , r^- Hence the modulus of (3.4), the m*^ term of the 
Dyson series, may be bounded above by 



\Xm\ < /3™ ^ • ■ ■ ^ |X^(rg,r2, . . . ,r^)| < trace 



»"i=l rm=l 



ml 



Noting that the zeroth term of the Dyson series is 
Xo = trace n,,v [Uge~^"] = trace n 



■fj^g-/3(H|?'-MJVv)' 



16 



we may re-sum the series to obtain 



trace -H, ,, [f/^e 



] - trace n^, 



< trace 



trace 



-^+(Wv) 



-/3g(l-M) 



m=l 



Thus 



V 



trace [f/^e-'^^^^-'^^^)] - trace [f/^e-^^-f^' 



(9) 



< 



e/39M(l - e-/^9) trace ^^(^^) [e"^^'] 



Since if,/ — /lA^'v — = \ J2x=i('^xi^x ~ 1) + 2A^^?T,a.) > 0, the second fraction is not greater 
than 1, implying 



Kiq)-c^viq)\< 



e/5</M(l -e" 



V 



which goes to zero in the limit V —> 00, as desired. 



□ 



3.3 Proof of Lemma 3.3 



The modified cycle density c^(g) may be simplified as follows: 

V V 00 



-/9(i-M)g 



^■■■^^{ii,i2,---,iq]k\ 



il = l iq = l k=l 



exp I -/3 (^dT\hy) + A J2{n, + N,){n, + iV, - 1) - /i ^ j | 



-I3(l-n)q V V 00 



V 



l^J^ X] ■ ■ ■ X] • • • , iq\Uq\il, 22, • • • , iq) 



i\=l in=l k=l 



{k\ exp -/3 j dT{hv) + \^^{n.x + + ^^^ij,x - 1) - 

x=l j=l j=l 



\k) 



x=l 
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and as (zi, «2, • • • , iq\i2, ■ ■ ■ , iq, h) 7^ if and only if ii = 12 



iq '■= i then 



j=i fe=i 



exp ^(ir(/iv) + A^(n^ + qSi^){n^ + q6i^ - 1) -/i^n^^ | 



I A:) 



^ 2^ trace 



i=l 



V 



V 



exp< -/3 dT{hy) + X{ni + q){ni + q - 1) + A^na.(na; - 1) -//^ 



x=l 



trace 



exp <( - /3 dVQiy) + A(ni + g)(ni + g - 1) 



+ A ^ n^{n^ - 1) - /i n^i. 



a;=2 



x=l 



since the trace is independent of the basis chosen. Hence we obtain 

^trace j-^CH^) exp { - (3{2\qni + Hy - yuiVy)} 



trace exp { - P{Hy - fiNy)} 



4 Proof of Theorem 2.3 



□ 



For convenience, denote Hg y := 2Xqni + (g — fi)q + Hy — fiNy. Due to Lemmas 3.2 and 3.3, 
we need to consider 

trace j-+(^^)[e-^-f^'>^-] 



c^{q) = hm 



v^oc trace j-+(K^)[e-'^^o.^] 



Recall that 



1 ^ ^ 



x,y=l 



x=l 



Motivated by the occurrence of the site-specific operator ni in the numerator of the expression 
of Cy{q), we expand the expression for Hy — fiNy to isolate the operators which apply to the 
site labelled 1: 

Hy - fiNy = (1 - fi)ni + Xni{ni - 1) + y - ^^al - y^a^ + Hy 

Xf^l x^l 



where 



Hy = {I - fi)^n^ - ^ Yl a>y + A5]]^xK-l) 



x=/=l 
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is an infinite-range-hopping Bose-Hubbard Hamiltonian for Av \ {1}. By denoting 
hq^y = (1 - + q) + A(ni + q){ni + q - I) + f 

we may tlien write 

Note tliat hq y — hq on J-'+(C) as V — )■ oo. 

We intend to completely segregate the Hamiltonian Hq y into two individual parts, one which 
operates solely upon the site labelled 1, and the other which applies only to the remaining 
V — 1 sites. What prevents us from doing this immediately is of course the "cross-term" 

— y^a* + — y^a,,. 

Y ^ y ^ 

x^l x^l 

Motivated by the Approximating Hamiltonian technique, we shall substitute this term with 

for a certain co-number R. Without loss of generality, we may take i? to be a non-negative 
real number. Fixing 

hq,v{R) = hq^v - R{ai + a*) 
then the resulting newly approximated Hamiltonian may be expressed as 

H^;:;{R) = hq,yiR) + Hy. 

In the arguments that follow, we shall either take R to equal in the variational principle, 
or a variable depending on V which tends to in the limit. 



4.1 Case 1: values of ^ such that = —he absence of condensation 

First we shall state and prove the following: 

Proposition 4.1 For all fi &W such that = 0, 

j.^ trace exp { - l3H^^^''{0)} _ ^ 

v^co trace j-^fH^) exp { - (3Hq^y] 

Proof: Using the Bogoliubov inequality: 

{A-B)b <ln trace - In trace < {A - B)a (4.1) 

for any R we obtain 



< In trace exp | — [3Hqy^ — In trace exp | — I3H^^ [R)^ 



</3(a,(±^^-i?) ) +h.c. (4.2) 
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Taking the left-hand side, since H^^{R) is a sum of two Hamiltonians which act upon 
different Hilbert spaces, traces and therefore expectations may be easily de-coupled, so one 
may see that 



V 



hq,v(R) 



V 



0, 



which is zero since Hy is a gauge-invariant Hamiltonian: ( Xlx^i ^*x) ~ ~ ^■ 



Hv 



On the right-hand side, note that (oi)// y = 0, also due to gauge invariance. Therefore we 
may simplify (4.2) to obtain 



hq,v{R) 



< In trace exp { — (3Hq^v} — In trace exp { — (3Hg^{R)^ 



</3(ai 



V 



+ h.c. (4.3) 



H. 



q.v 



For this case we shall take R = = 0. We therefore obtain 

< In trace exp { - f3Hg^y} - In trace exp { - /3H^^y''{0)} < (3 (a^ 

Now by the Schwarz inequality 



V 



H. 



+ h.c. 
(4.4) 



Xlx^l ^*x 

V 



H. 



T.x< 

V 



H. 



q,v 



q,v 



V 



To consider this let H^ y := Hq^y + s'^'^^^^^^ and H^y := H^y — 2Xqni. Using the Bogohub 
inequality (4.1), with A = —/3H^y and B = —l3H^y, so that A — B = 2/3Xqni, one has 



ov 



2f3Xq{ni)H^^^ < In trace e-^-^^V - In trace e'^^^^V < 2/3Ag(r2i)|j. ^ = 2/3Ag 



^q.v 



q,v 



V 



(4.5) 



The last equality is due to the fact that the system H^y is invariant under permutation of 
the sites of the lattice. This identity implies the following: 

(i) With s = we get (rii)// ^ < {N)Hy/V, thus {ni)H y is bounded and in the limit 



{ni)i 



V 



^ ^ 0. 



(ii) 



< In trace e ^^^.v _ In trace e ^^^V < y y 



2/3Xq (^)Hiy 
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indicating that in the hmit, the pressures are the same for the two Hamiltonians. 
By using Griffith's Lemma we see that the condensate densities (the derivatives with 
respect to s at zero) are both equal to zero (since we are considering the case = 
here). That is 

Using these facts, one sees that the right-hand side of (4.4) goes to zero in the hmit and we 
can conclude that 

tracej- (^^)exp{ -/3ifq,v-} 

lim 7-^ T^T^ = 1. 

v^^ trace exp { - PH^^/{0)} 

□ 



Using Proposition 4.1, one immediately obtains the desired result: 

<="(,) = lim <,(,) = lim *--'"'^.m., "-A 

v^oo v^oo trace t+chv) exp { - pH^^y \ 

trace^^(^,,)exp{-/3i7^^7(0)} 

= lim ) — —4- 

v^^ trace ^^(^^) exp { - /3//oT(0)} 

- lim ^^^^^^-^+(^^-1) i ~ f^Hv] trace j-(c) exp { - I3hq^v] 
trace j-+CH^,_o exp { - (3Hy] trace j-(c) exp { - /3/io,v} 

_ trace j-(c) exp { - /3/ig(0)} 
trace ;r(C) exp { - /3/io(0)} 

in the limit since hq^v = ^,v(0) — /ig(0) on J^(C). 



4.2 Case 2: For any /i 

Considering the case when > 0, if we insert i? = in the constraint inequality (4.3) then 
its left-hand term is strictly negative, but its right-hand term is strictly positive, rendering 
the previous argument useless here. 

We therefore introduce a gauge-breaking term v + ^ '^*x i^^o the Hamiltonians Hq^y 

and H^y^{r^). Without loss of generality we may assume u to be real and positive, so denote 

V 

Hq,y{v) = Hq^y - V ^(o^ + <) 
x=\ 

and its corresponding approximation as 

V 

E^y\R, v) = H-J{R) - u Y.(a^ + <)■ 

x=l 
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Again we wish to separate this Hamiltonian into parts, one acting upon the site labelled 1, 
the other on the remaining V — 1 sites. If we define: 

hq,v{r, y) := hq^v{r) - v{ai + a\) 

= (1 - + g) + AK + + g-l) + ^-(r + v){ai + al). 

and 

Hy{v) = Hy - + al) 

then we may write H^J^{R, u) = hq v{R, i^) + Hv{u). Denote limy_i.oo hg v{r, u) = hg{r, u) on 
J-'+(C), i.e. obtain a "gauge symmetry broken" single site Hamiltonian 

hq{r, v) = {I - fi){n + q) + X{n + q){n + q - 1) - {r + iy){a + a*). 

We shall first prove the following proposition: 



Proposition 4.2 For each v > 0, there exists a sequence {uy E M : Uy E [u, u + 1/v^]} 
independent of q such that 

lim t'^acej-+(^^)exp{-/3iJg,^(z/v.)} _^ 



v^oo trace jr^(^uv) exp {-/3H^pp{r^{uy), Uy)} 

where r^{v) > zs the non-zero solution of 

2r = (a + a*)/,^(,,^). (4.6) 

Note that \im.y^Qr^{v) = r^, the maximal solution of {2.5), i.e. the positive square root of 
the condensate density. 



Proof: There is no immediate correlation between the chosen R and as yet. For each 
z/ > take a sequence Vy which tends to as V — oo. Then using the Bogoliubov inequality 
again, we obtain 

\ V y J I H^l^{R^Vy) 

< In trace exp {—f3Hq^y{h'y)} — In trace exp |— i/^)} 

<f3(aJ?^-R]) +h.c. (4.7) 



V 



As above, the left-hand side may be reduced to 



Ex^i p^ \ - In \ I ^^^^ p 



^ y / K^P{R,Uy) \ ^ ' Hy(vy) 



If we replace R with the term 



^ ' Hy{uy) 
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then the left-most side of (4.7) is identically zero and we get that 



< In trace exp {—(3Hq^v{i^v)} — In trace exp [—(3H^^^ 



Sf^,v{'^v),'^v 



)}• 



Hence 



lim inf 



trace exp {—/SHq^yluy)} 



v^oo trace exp { — 



)} 



> 1. 



(4.9) 



Now considering the right-hand side of (4.7) for any fi. Using the Schwarz inequality as 
before: 



<{ni) 



a. I - R 



V 



Hq,vM \\ y 



1 



1 

R \ ( ^^-R]y + 



V 



{ni)H,y{v) + 1 
V 



where we have taken 



1/2 



V 



(4.10) 



Again in order to consider this, insert Hq^v{v) and Hqv{v) := Hq^y{v) — 2\qni into the 
Bogoliubov inequality, to obtain 



2l3\q{ni)H^^y[y) < In trace e-'^-^''^^'')- In trace e-^^«-^(^) < 2(3Xq{ni)^ = 2/3 Ag 



which implies the following facts: 



V 

(4.11; 



(i) {ni)H, ^ H,y-{v)l^ ^ and hence 



(ni> 



gq,v('^) 



V 



— )■ as y — )■ 00. 



(ii) Since Hq^y{u) = (g - /i)g + Hy{v) - fiNy, 
< 



/3(g - /i)g ^ 1 1, trace e-'^^^-^^)-^-) - 1 Intracee'^^-H < '^^^ ^^^^^^^^ 



V 



V 



V 



V V 



with which one may show that in the limit, the pressures are the same for Hy{v) and 

hm Py[Hy{v)] = hm Py[Hq^y{v)]. (4.12) 



Using fact (i) from above, we find that the only term we need yet be concerned on the right 
hand side of (4.7) is the first term of (4.10), whose behaviour in the thermodynamic limit is 
still unknown: ^ 
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To deal with this we shall take R to be the following: 

<vM = ^(E«-/ (4.13) 

SO that one has 

^o = 47fE«^-(E«^) V (4.14) 
Now we shall state and use some lemmas, which are proved later: 

Lemma 4.1 For fixed z/ > 0, a positive integer q, and (5o is defined as (4.14), then there 
exists a sequence {uy G M : i^v- ^ K + independent of q which tends to v as 
y — )■ oo, such that for large V we have the approximation 



w 



for some constants u, w and B, independent of u and q. 

Using this lemma, for large V and fixed q, we obtain the following estimate 

const) 

In trace exp {-^Hg^^^Uy)} - In trace exp {-f^H^^/^r-^^^^Uy), Uy)} < -y^j^ 

where Uy ^ u as V ^ oo, implying that 

trace exp {-/3i7g_,,(i/v-)} , , 

limsup p , — TY < 1- (4-15) 

v^oo trace exp [-l3H^^p{r'^^y{i'y),Uy)\ 

Lemma 4.2 For a fixed u > and for any sequence {i^y} which tends to u as V ^ oo, then 

(uy) = lim r+ (uy) = r^(z/) 

where r^{v) is the unique non-zero solution to the Euler- Lagrange equation {4-6). 

For clarity, it is best to use the following short-hand for this argument: 

ay = trace jr^-Hv) exp { - (3Hg^y{uy)} 
by = trace jr^-Hy) exp { - PH^J^r^^Uy), Vy)] 
Cy = trace jr^ny) exp { - l3Hg^/{r~^{iyy), Uy)] 
dy = trace jr^ny) exp { - /3iJg';7(^M,v('^v), z^v)} 



The penultimate step is to prove the following: 

lim — = 1 and lim = 1. 

y-S>00 Cy V^OO dy 
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Considering the first, note that 

by _ trace ^(oe-'^'^^^'^''^'^^)'''^) 
Cv ~ trace ^(c)e-^'''(Vv(-v),-v)' 

Once again by the Bogohubov inequahty (4.1), with A = —f3hg{r~y{h'y),h'y) and B = 
— Phqir ^{uy) , Uy) then A — B = f3{r^{iyy) — 'r^^y(z^v))(a + a*) and we obtain 

(^{r^.M - r-,,(i^v)) (a + (^*) h,[r ,(yy\vy) 

<\Tyhy- Inc, < /3(r^(z/v) - r-,(z/,)) (a + «*)^(.-^(.,),.^) • (4.16) 

By continuity one may see that hmy_^oo '''/^(^^v) = "^iiiy)- Then as V — oo, both the left and 
right hand sides of (4.16) go to zero, implying the first result. A similar procedure may be 
used to show the second. 

We want to show that limv-^.oo = 1- Since (4.15) implies that limsup ^ < 1, we have: 

a.r ^ lim supv^nr, §^ 
limsup — = limsup -r^ < r — — < 1. 

V^oo by V^oo limy^oorf^ 

The infimum limit follows similarly from (4.9), proving the proposition. □ 



With the assistance of Proposition 4.2 we then have our result: 

lim <(g,..) = hm trace ^,(^,) exp {- /3g„.(..)} 
v^oo v^oo trace r+{Hy) exp { - (3HQ^y{uy) \ 

_ ^.^ trace j-^cHy) exp { - l3H^^/{r^^y{uy), Uy)} 
v^oo trace jr^(^uy) exp { - l3H^^/{r^^y{iyy), Uy)} 

_ j.^ trace j-^cHv,_o exp { - iSHyjuy)} trace j-(q exp { - l3hq^y{r^^y{uy), Uy)} 
trace j-_^(^^„_^) exp { - PHy{iyy)} trace j-(c) exp { - l3hQ^y{r^^y{iyy), Uy)} 

_ trace j-(c) exp { - f3[hg{r^{u)) - u{a + a*)]} 
trace j-(c) exp { - /3[/io(r^(i^)) -y{a + a*)]] 



4.3 Proof of Lemma 4.1 

Proof: In the Appendix we prove that (see (A. 6)) there exist constants u and w such that 
for fixed G M, all v < and all g G N, 

H,.yi.) < UA^) + + ^ (4.17) 

with fq,v{y) '■= {So,So)Hgy{u) where (■,■)// is the Duhamel inner product, see (A. 5). Set 
Co = ^ E «x. 

One may check that for > 0: 
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This follows from the fact that when e C, one can write fq,v{i^) = dyd^pvlHg^yli/)], but 
since pvlHg^Yi^)] does not depend on the argument of p we can use polar coordinates to get 
(4.18). We want to show that limv_s.oo fq,v{^)/V ~^ 0. Consider (4.18), multiply both sides 
by u and integrate: 



for [v, V + \/\/V] C [0, uq\. Using the bound (A. 8) proved in the Appendix, then there is a 
constant B such that for all g G N we obtain 

j v'fq,vW)du' <b[2u + -j=j . (4.19) 

Now let 

oo 

F,.(i.) = J]e-%,v-(^). (4.20) 

9=1 

In the Appendix (see (A.7)-(A.9)) we show that there exist constants a and h independent 
of q such that for all v < vq and g G N: 

fqA^)<- + bV. 

V 

Therefore the series (4.20) is uniformly convergent in v and since each term is continuous in 
Fy(v) is also continuous. From (4.19) we obtain: 



By the Mean- Value theorem, there exists some Vy G [i^, + 1/ \/V\ (independent of q) such 
that 

v'Fy{v')dv' = -^Fy{Uy) (4.21) 



which implies that 

Fy{Uy) < B (2VV— + — 

For any positive integer g, since e~'^fq,v{iy) < Fy{u), then 



fq,vM < Be" (2VV— + —] . 



Thus we have a sequence {uy ^"E. : Uy ^ [u^u ^- 1/yV]} satisfying (4.21), independent of g, 
which tends to z/ > as — oo, such that for large V we have the estimate: 

Combine this with (4.17) to complete the proof. □ 
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4.4 Proof of Lemma 4.2 

Proof: Before proceeding, we need to show the following: for fixed z/ > 

lim pv[Hy{u)] = hm pv[Hv{iy)], 

i.e. a single site's contribution is irrelevant in the thermodynamic limit. Recall that 

The corresponding pressure may be expressed as 
Vv^v\ = -^lntracej-_^(^^)exp 

1 r / 1 

= ^ In trace exp | - ( ^ a^) (a^ - + ""^("-^ ~ 



+ (?(FrT)-''jg"- 

\ 2 

^" Pv-i[i^v-i(/3(V),A(^),/x 



V-1 



y I Fv-Li^^v-i\y \ Y )-,'^\v-l)'t^ V(V-1) 

where abusing notation temporarily we have explicitly included the parameters of the IRH 
Bose-Hubbard Hamiltonian, i.e. we write the pressure of (2.2) as pv[Hv\ = Pv[Hv{l3, A,/i)]. 
Then in the limit, with the use of the Bogoliubov inequality, one may verify that 



lim^p^Hy] = ^im ^ j [H^.ii(3 (^) , A (^) ,/i - ^(^)] 



lim Pv[Hv] =p{l3,^). 



Now proceeding to prove this lemma, recall that we chose 



\ y 



where Hy{v) = Hy — v Xl^^il'^a: + ^*x) is the gauge-broken IRH Bose-Hubbard Hamiltonian 
on all sites of the lattice barring the site a; = 1. Fixing a value of z/ > 0, there exists a unique 
r^jyv) > as the solution to the Euler-Lagrange equation (4.6), i.e. 

'^r^iv) = (a + a*)^(rvH,i.). 

The pressure pv[Hv{i')] is convex in u and its thermodynamic limit is differentiable for all 
u > 0. By Griffith's Lemma, we have 

lim ^p,[Hy{uy)] = ^ hm py[Hy{u)]. (4.22) 
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The left hand side of this evaluates to 

d ~ 1 2 

As shown above, we have that 
lim py[Hy{u)] = lim pv[Hy{i')] 

V^oo 1/— s>oo 

= -r^{uf + i In trace exp - l)n - \n{n - 1) + + iy){a + a*)] }. (4.23) 

so the right-hand side of (4.22) will become (also using (4.23)) 



lim pv[Hv{u)] 
d 

du 



"^iii^Y + ~B In trace exp (/3[(/i — l)n — Xn{n — 1) + {r^{h') + v){a + a*)]) 

, ^dr,,(u) f dr,,(u) \ , 
= -'^'■M^ + + Ij (a + a*),^(.^M,.) 

= ((a + - 2r^(i/)) + (a + a*)/,^(^^(^),^) 



as desired. 
Similarly, taking 



where Hq^^{y) = Hg y — v Y^x=\^^x + ^*x)- Label the corresponding pressure for this Hamil- 
tonian as pv\H q^,r(y^. Recall the expression (4.12) that we previously derived: 

lim py{Eq^y{y)\ = lim py[Hy{u)]. 
As above by Griffith's Lemma, we have 

lim r+(z/v.) = ^ lim -J^pylHg^yiuy)] = Jim Pv[Hy{iy)] = r^(z/) 
as above. □ 



5 Proof of Theorem 2.4 



For those values of /i such that = 0, i.e. in the absence of condensation, first note that 
the density (from (2.4)) may be expressed as: 

d d I 

" " 5/i ^^^^^^ ^^^^ ^ ~ ^^^^ ~ ^ ^^^^ ~ ^ 

trace jr(-c) [nexp | — /3[(1 — fi)n + \n{n — 1)]}] 
trace exp | — /3[(1 — ji)n + \n{n — 1)]} 
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Similarly, from Theorem 2.3 one immediately obtains: 

trace j-(c) exp { - - + g) + \{n + q){n + q - 1)]] 



trace jr(c) exp { — — fj.)n + \n{n — 1)]} 



Label the denominator H := trace jr(c) exp{— — /x)n + \n{n — 1))}. The operator n in 
this context counts the number of particles on the site, so in terms of a basis of occupation 
numbers, it has eigenvalues /c = 0,l,2,.... Summing over this basis 



^ _j oo oo 

c^{q, 0)=^J2Y. { - ^[^(1 -f^)ik + Q) + Hk + q){k + q - I)]} 



q=l q=l k=0 

and shifting the sum 



oo oo 



q=l k=q 
^ oo k 

=^ ^ ^ exp { - m - f^)k + Xk{k - 1)] } 

k=l q=l 
^ oo 

=^ k exp { - - fi)k + Xk{k - 1)] } 



k=l 



=— trace j-(c:) [nexp { — — ii)n + Xn{n — 1)]}] 



Therefore the absence of condensation implies that the sum of all finitely long cycle densities 
equals the system density. □ 



Appendix A. Some useful inequalities and bounds 

By the operator inequalities CqCq < A''^, Ylix'^'x — "y^ 

uVv{co + c*) < u^clco + V < u^Ny + V, 

it is clear that the Hamiltonian with sources, Hq^vi^)-, is superstable for fixed g G N and 
A > 0, i.e. 

HqA^) > A-^ - (A + /i + z/2)iV, - V. 
Using the Bogoliubov inequality (4.1), since Hy^P, /i — 2gA, u) — Hq^v{(3, p, i^) — q{q — p) = 

2qX X]j=2 o'^^ ^^^^ 

i=2 I n^y{^,ii.,v) \ j=2 / HvW,fi-2q\,u) 
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Then we may find an upper bound (which is independent of q) for the expectation of the 
number operator with respect to Hq Y{v), using (4.11) and (A.l), as follows: 

Now 

y \ 

since {j^v) ^ (fSfiu) monotonically increasing in fx. Hence 

(Nv)^ r •.<(—) +(Nv)^ ,<2(Ny)^,,. (A.2) 



By the superstability of the Hamiltonian Hviy), there exists a function M > such that 

{^v) Hv(u)/'^ — ^ ^ ^'^d 1/ < i/q- Using the fact that {a^)HgvM ~ {^x)HgvM 

(A.2) we have 

for all u < Uq and g G M. 

Considering the term {Sq6q)jj ^^^^ stated in Lemma 4.1, we follow the procedure in Appendix 
1 of Bru and Dor las [7] to write 



where the Duhamel inner product ( ■ , ■ ) is defined as follows: 

r./3 



{A,B)h = -3^ / trace [A*e-^^-'^^ Bq-'^] ds 
Jo 



(A.5) 



with Z = trace e . Since Sq = Cq + (cq) one may then evaluate that 

[S;, [H,A^), So]] = K [H,A^), Co]] = -/i + A + 2A^^^^ 
to give that (A. 4) has an upper bound of the form: 

(<^o^o)//,,vH ^ (^0, So)H,^yiu) +u^ + w (A.6) 
(using (A. 3)) for some constants u and w (independent of q). 

We set fq,v{ty) = {So,6o)h vM- wish to show that fq,v{^) is bounded by a constant 
independent of q. For fixed V and u < uq, we have 

4/3/,,v(^^) = l-^pAHqA'^)] + ^PvlHqA'^)]- (A.7) 
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By (A. 3), we have for all g G N that 

d 



The second derivative term of (A. 7) is: 



^2 2 



using the fact that {A, A) < ^{A*A + AA*) 

< ((Co + C*)(Co + C*))^^ ^(^^ 



<V{AM+l + 2^/2Mi2MTT)) (A.9) 
by the Schwarz inequality, (A. 2) and (A. 3). 
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